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• j—j Abstract 

Phase contrast techniques, such as differential interference contrast (DIC) mi- 
croscopy, are widely used to provide morphological images of unstained biological 
samples. The trans-illumination geometry required for these techniques has restricted 
their application to thin samples. We introduce oblique back-illumination microscopy 
(OBM), a method of collecting en face phase gradient images of thick scattering sam- 
Qh pies, enabling near video-rate in vivo phase imaging with a miniaturized probe suitable 
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for endoscopy. 
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Phase contrast microscopy techniques are widely used in biological research because they 
can provide high resolution images of unlabeled samples even when these are nearly trans- 
parent. For example, differential interference contrast (DIC), which reveals lateral phase 
gradients, is one of the more popular techniques because it provides apparent 3D sample 
relief using a standard microscope equipped with a lamp and a camera fil. Even simpler 
techniques based on oblique illumination can also be used with standard microscopes and 
provide similar imaging as DIC [2]-[5] . However, in order to reveal lateral phase gradients, 
all these techniques must be operated in trans-illumination configurations, limiting their use 
to thin samples. Because one is often constrained to working instead with thick samples 
(e.g. endoscopy or in vivo applications), there is a clear need for a method that can provide 
DIC-like imaging in thick samples. We introduce such a method here. 

Our technique is called oblique back-illumination microscopy (OBM). As its name sug- 
gests, this technique is based on a similar principle as oblique illumination microscopy, with 
the notable difference that illumination and detection both occur from the same side of the 
sample (i.e. in a reflection geometry), allowing its application to samples of arbitrary thick- 
ness. Unlabeled samples can, of course, be imaged with microscopes based on direct light 
reflection. The most successful of these for tissue imaging is optical coherence tomography 
(OCT), which, like OBM, can also be operated in a widefield en face configuration rcjj. How- 
ever, microscopes based on light reflection intrinsically reveal only sample structure that 
varies rapidly in the axial direction, such as sharp interfaces or particles much smaller than 
the illumination wavelength Yn. In contrast, microscopes based on light transmission are 
not subject to this constraint, and can reveal even slowly varying sample structure in the 
lateral direction, thus providing images of subtle sample features impossible to see in reflec- 
tion mode. An important characteristic of OBM is that, even though it is configured in a 
reflection geometry, it is, in fact, a transmission microscope in disguise. In effect, OBM uses 
multiple scattering in tissue to convert epi-illumination into trans-illumination. Because the 
illumination source is offset from the detection optical axis, the trans-illumination is oblique 
(Fig. [I]). While illumination obliquity directly leads to phase gradient contrast, image in- 
tensity is also influenced by sample absorption. The use of two off-axis sources diametrically 
opposed to one another permits the acquisition of two raw images with similar absorption 
contrast but with phase gradient contrast of opposing sign (Supplementary Fig. [I]). The 
subtraction of these raw images enhances phase gradient contrast while canceling absorp- 
tion contrast; addition of the raw images has the opposite effect, revealing only absorption 
contrast while canceling phase gradient contrast. By this method, the sequential acquisition 
of two raw images using alternating illumination sources decouples absorption and phase 
gradient contrast. 

We present results obtained from a miniaturized OBM built using a flexible endomi- 
croscope probe comprising a distal micro-objective and an imaging fiber bundle comprising 
30,000 fiber cores (Online Methods). Illumination from two independently controlled light 
emitting diodes (LEDs) was delivered to the sample via two optical fibers attached opposite 
one another to the micro-objective housing. The endomicroscope probe was designed to 
operate in contact mode, meaning that light reflected directly from the sample surface is 
not collected by the micro-objective. Instead, the micro-objective only collects light that 
has been multiply scattered in the sample and redirected upwards through the focal plane, 
located here at a depth of 60 /im (the working distance of the micro-objective). Unless 
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stated otherwise, all images presented here are individual frames from movies acquired and 
displayed at a net rate of 17.5 Hz using a double-shutter camera that reads images pairwise. 
The exposure time per raw image was typically 1 to 5 ms. 

Figure [2] shows a phase gradient image of a 45 fim polystyrene bead embedded in a 
scattering tissue phantom consisting of 2 fim beads mixed in agarose. The phase gradient 
induced by the bead was observed to be approximately linear. Note that since the phase 
gradient image is derived from a difference of raw images, it contains both positive and 
negative values (zero is represented by an intermediate gray level) . Because phase gradients 
must, by definition, arise from apparent sample structure, they must also arise from the 
vicinity of the focal plane (objects out of focus are blurred and exhibit little structure). 
Thus, phase gradient imaging exhibits apparent out-of-focus background rejection. This is 
manifest in Supplementary Video [TJ in which only the in-focus 2 fim beads are visible. 

Figure [2] also shows simultaneously acquired absorption and phase gradient images of 
the chorioallantoic membrane (CAM) in a day 11 chick embryo in vivo and in ovo. Note 
the low contrast of the absorption images compared to the phase gradient images. Because 
the absorption images are derived from the sum of raw images, their values are positive 
definite (zero is black), and they do not exhibit out-of-focus background rejection (such 
imaging is similar to orthogonal polarization spectral imaging (81). Compared to images 
obtained from OCT, OBM images are speckle-free. Moreover, owing to the large photon 
fluxes involved, they are also relatively shot-noise free. Intensity noise arises, in our case, 
mostly from inhomogeneous image sampling due to an uneven distribution of fiber cores in 
the imaging fiber bundle (Online Methods). 

Finally, Figure [3] shows simultaneously acquired absorption and phase gradient images 
of excised, unstained mouse intestinal epithelium. The absorption images are essentially 
featureless in this case, whereas the phase gradient images are information rich. For example, 
crypts of Lieberkiihn and ileal villi are readily visible in Figures [3^ (Supplementary 
Video [6]) and [3ji (Supplementary Video [7]), respectively, illustrating the potential of 
OBM for in situ histopathology and "optical biopsies" . 

In summary, we have presented an apparatus that, to our knowledge, is the first to provide 
sub-surface, DIC-like (i.e. transmission-like) phase gradient imaging from thick scattering 
tissue in a reflection geometry. The apparatus is simple, fast, robust, and inexpensive, 
making it broadly appealing to biological and clinical researchers alike. 
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ONLINE METHODS 

Hardware setup 

White light from two LEDs (Luxeon Star MR-WC310-20s) was coupled into optical fibers 
(Thorlabs BFL48-1000; 0.48 NA; 1000 /xm core) using aspheric condenser lenses (Thorlabs 
ACL5040-A). Illumination light was launched by the fibers into the sample (25 mW per 
channel at the fiber output), where it was redirected through the focal plane by multiple 
scattering and collected by a micro-objective (Mauna-Kea Technologies; 2.6 mm diameter; 
lx or 2.5 x magnification; 60 /im working distance; water-immersion; 0.8 NA) coupled to a 
coherent imaging fiber bundle (30,000 cores; 600 /xm active area). The separation distance 
between the fiber and the micro-objective probe was approximately 1.8 mm. The proximal 
face of the fiber bundle was imaged with standard microscope optics (Olympus Plan 10 x 
0.48 NA air objective, Linos AC / = 200 mm tube lens; 4/ configuration) and recorded with 
a digital camera (PCO Pixelfly USB; 14-bit; 2x2 binning; 35 fps; 1 — 5 ms exposure time 
per illumination direction) . The camera was operated in double shutter mode to reduce the 
inter- frame delay between exposures (200 /xs), minimizing motion artifacts (9). Illumination 
power delivered by the left and right optical fibers was triggered (Thorlabs LEDD1B) to 
overlap with the first and second frame in the each image pair, respectively. Frame rate 
was limited by the camera readout time. Image acquisition and display was performed 
using custom written software (National Instruments Lab VIEW 11.0). Illumination gating 
and camera exposure were synchronously controlled using a data acquisition card (National 
Instruments PCI-6221). 

Image processing 

A preprocessing routine described previously |9] was first used to correct for the quasi- 
periodic sampling pattern imparted by the fiber bundle cores. Each raw image was then nor- 
malized by its respective low-pass filtered version (Gaussian filter kernel with a = 80 pixels) 
to correct for non-uniform illumination profiles and thus "flatten" the images. The two nor- 
malized images were then either added or subtracted to produce absorption-only or phase 
gradient-only images, respectively (Supplementary Fig. [I]). Image processing was per- 
formed with a graphics processing unit (NVIDIA GTX280) using custom-written software 
written in CUDA-C [To] . 

Monte Carlo simulations 

CUDAMCML [II], a modification of MCML [12] enabling execution on graphics processing 
units (GPUs), was used to perform the simulations. CUDAMCML was further modified 
to execute on a cluster of CUDA-enabled workstations [13]. A semi-infinite slab geometry 
was modeled with tissue optical parameters n t i ssue = 1.37, l s = 150 /im, Z* = 3000 /xm and 
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g = 0.95 (n is index of refraction, l s and I* are the scattering and transport mean free 
path lengths, respectively, and g = 1 — l s /l* is the anisotropy factor). Illumination fiber 
parameters were n^er = 1-37, diameter = 1000 /im and numerical aperture, NA = 0.48. 
Micro-objective probe parameters were n pro b e = 1.37, diameter = 240 /xm and NA = 0.8. 
Fiber-probe separation was d = 1818 /im. A Henyey-Greenstein phase function was used 



to characterize photon scattering events 14 . 10 photons were processed to estimate the 



distribution of exit angles of the detected photons as a function of fiber-probe separation 
(Supplementary Fig. [2]). Both the total detected intensity and median exit angle were 
observed to decrease with increasing fiber-probe separation. 10 5 photons were processed 
to estimate photon path density as a function of lateral position and depth, revealing the 
so-called photon banana (Fig. ^>). 



Tissue phantom preparation 

The scattering tissue phantom was prepared by heating a 30 mL solution of 2% (w/v) agarose 
(Sigma A5093-100G), 5%2/xm diameter polystyrene beads (Polysciences 19814-15), and 
0.1%45/xm diameter polystyrene beads (Polysciences 07314-5) in H 2 to 75°C on a hot- 
plate, followed by pouring the mixture into a 60 mm x 15 mm cell culture dish (Corning 
430166). The phantom was covered with paraffin film and left to cool to room temperature 
before imaging. The optical properties of the bulk medium were l s = 74 /im, /* = 1040 /xm 
and g = 0.93, as estimated using Mie theory. The indices of refraction of hydrated agarose 
gel and polystyrene beads were n = 1.35 and n = 1.59, respectively |15|. Imaging was 
performed through water. 



Chick embryo preparation 

Fertilized Gallus gallus eggs (Carolina Biological Supply Co. 139290) were stored in an in- 
cubator at 37°C and 50% humidity, being turned every 7h to prevent fusion of the chorioal- 
lantoic membrane (CAM) with the shell membrane. Imaging was performed at embryonic 
day 11. A 1 cm diameter region of the shell and shell membrane was removed exposing the 
embryo and CAM. A layer of 37° C saline was dripped over the preparation before imaging in 
ovo with the OBM probe. Following imaging, the embryos were euthanized by hypothermia 
by storing the eggs at — 15°C. 



Mouse tissue preparation 

Six week old C57 black 6 mice were euthanized by CO2 inhalation and the gastrointestinal 
tract was immediately excised and washed with 4% paraformaldehyde. The colon and small 
intestine were cut longitudinally, unrolled and cleared of fecal matter. The preparations were 
stored in 4% paraformaldehyde for several days. Before ex vivo imaging, the tissues were 
pinned to a silicone elastomer slab (Sylgard® 184, Corning) to expose the apical surfaces. 
Residual mucus and fecal matter were gently washed away with saline before imaging. The 
animals used in this study were treated in accordance with the guidelines of the Institutional 
Animal Care and Use Committee of Boston University. 
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Lateral position (jam) 

Figure 1: An OBM setup with a contact-mode endomicroscope probe, (a) Illumination from 
two LEDs is sequentially projected into a thick sample by diametrically opposed optical 
fibers attached to the probe housing. Multiple scattering in the sample redirects the light so 
that it trans-illuminates the focal plane of the probe micro-objective (magenta dashed line). 
An image from the focal plane is then relayed by a flexible fiber bundle and projected onto 
a digital camera, (b) Close-up of the probe distal end, onto which is superposed a density 
map obtained by Monte Carlo simulation of the light energy in the sample that was injected 
by a single fiber (from left) and collected by the micro-objective (Online Methods). Note the 
obliqueness of the light distribution through the focal plane, (c) Oblique trans-illumination 
is partially blocked by the micro-objective back aperture. Index of refraction variations at 
the focal plane refract the light causing changes in the image intensity that are proportional 
to the slope of the variations, thus leading to phase gradient contrast. 
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Figure 2: Demonstration of OBM in a tissue phantom and in vivo, (a) Phase gradient image 
of a 45 fim polystyrene bead in scattering medium (2 fim beads in agarose, Supplementary 
Video Note 2 /im beads are readily visible, (b) Corresponding phase gradient profile, 
(c-f) Simultaneously acquired absorption (c,e) and phase gradient (d,f) images of the CAM 
vascular system of day 11 chick embryo. Individual red blood cells (RBCs) and vessel walls 
are clearly visible (Supplementary Video [2]). (f) Moving RBCs are highlighted in red 
using a sliding 3-frame temporal variance filter (Supplementary Video [3]). (g) A CAM 
vasculature mosaic reconstructed from Supplementary Videos [4] and [5} Scale bars are 
(a, c-f) 20 /im, and (g) 50 /im. 
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Figure 3: Demonstration of OBM in excised mouse intestinal epithelium, (a-f) Simultane- 
ously acquired amplitude (a,c,e) and phase gradient (b,d,f) images taken with a lx micro- 
objective (field of view (FOV) 600 /mi). (g,h) Higher magnification phase gradient images 
of the epithelium in the distal colon (g) and small intestine (h) were taken with a 2.5 x 
micro-objective (FOV 240 /im). Crypts of Lieberkiihn (cr), crypt lumens (cl) and lamina 
propia (lp) are indicated with arrows (g, Supplementary Video [6]); ileal villi (v), colum- 
nar epithelium (ce) and lamina propia (lp) are indicated with arrows (h, Supplementary 
Video [7]). Scale bars are (a-f) 75 /im, and (g,h) 30 /im. 
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Supplementary Figure 1 


Comparison of added versus subtracted raw OBM images 


Supplementary Figure 2 


Photon exit angle distribution estimated with Monte Carlo simulation 


Supplementary Video 1 


Manual focusing through scattering tissue phantom 


Supplementary Video 2 


CAM vasculature and demonstration of axial resolution 


Supplementary Video 3 


CAM vasculature with moving RBCs highlighted in red 


Supplementary Video 4 


Comparison of absorption versus phase gradient images in CAM 


Supplementary Video 5 


Phase gradient mosaic of CAM vasculature 


Supplementary Video 6 


Morphological features of mouse distal colon 


Supplementary Video 7 


Morphological features of mouse small intestine 




Supplementary Figure 1: 45 jim polystyrene bead suspended in scattering tissue phantom. 
(a,b) raw images under oblique back-illumination from two opposing directions, (c) Addition 
of (a) and (b) cancels phase gradient contrast and emphasizes absorption, (d) Subtraction 
of (a) and (b) cancels absorption contrast and emphasizes phase gradients. The 2 /im beads 
used to build the tissue phantom are readily visible only when they are in focus. Scale bars 
20 /im. 



10 



Supplementary Video 1: OBM exhibits apparent axial resolution, as is demonstrated by 
focusing through suspended polystyrene beads. Scale bars 20 /im, imaging speed 5 Hz. 
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Supplementary Figure 2: Monte Carlo simulations estimate photon exit angle distributions 
at different fiber-detector separations. The exit angle corresponds to the tilt angle of the 
detected photon's path relative to the micro-objective optical axis (positive angles point 
away from the source). Five fiber-detector separations were considered: 1830, 1730 and 
1910 /xm correspond to the middle, left, and right extremes of the 2.5 x micro-objective 
FOV, respectively, while 910 and 3970 fim correspond to positions roughly half and twice 
these distances, respectively. The median exit angle for each distribution is noted with a 
dashed line. A Lambertian exit angle distribution, characteristic of isotropic illumination, is 
also shown for comparison. 
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Supplementary Video 2: CAM vasculature of day 11 chick embryo visualized with OBM. 
Capillary vessel walls are clearly visible, as are the dynamics of individual red blood cells. 
Axial resolution is made apparent by manually focusing up and down between the CAM 
mesoderm and ectoderm. Our imaging depth is limited here by the 60 /im working distance 
of our micro-objective. Deeper imaging could be achieved with a longer working distance. 
In principle, maximum imaging depth of OBM is expected to be similar to that obtained 
with DIC, namely on the order of l s . Scale bar 20 /im, imaging speed 17.5 Hz. 



Supplementary Video 3: OBM video of CAM vasculature of day 11 chick embryo. Moving 
RBCs are highlighted in red using a sliding 3-frame temporal variance filter. Scale bar 50 /im, 
imaging speed 17.5 Hz. 
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Supplementary Video 4: Simultaneously acquired absorption and phase gradient OBM video 
of CAM vasculature of day 11 chick embryo. The probe was scanned over the sample using 
manually controlled translation stages. Scale bars 20 /im, imaging speed 17.5 Hz. 




Supplementary Video 5: Mosaic created from OBM video of CAM vasculature of day 11 
chick embryo. Scale bar 50 /im, imaging speed 17.5 Hz. 
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Supplementary Video 6: OBM video of crypts of Lieberkiihn in excised mouse distal colon. 
Scale bar 30 /im, imaging speed 5 Hz. 




Supplementary Video 7: OBM video of ileal villi in excised mouse small intestine. Scale bar 
30 /im, imaging speed 5 Hz. 
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